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ABSTRACT 

It has been shown that the observed temporal distribution of transient events 
in the cosmos can be used to constrain their rate density. Here we show that 
the peak flux-observation time relation takes the form of a power law that is 
invariant to the luminosity distribution of the sources, and that the method 
can be greatly improved by invoking time reversal invariance and the temporal 
cosmological principle. We demonstrate how the method can be used to constrain 
distributions of transient events, by applying it to Swift gamma-ray burst data 
and show that the peak flux-observation time relation is in good agreement with 
recent estimates of source parameters. We additionally show that the intrinsic 
time dependence allows the method to be used as a predictive tool. Within the 
next year of Swift observation, we find a 50% chance of obtaining a peak flux 
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greater than that of GRB 060017 - the highest Swift peak flux to date - and 
the same probability of detecting a burst with peak flux > 100 photons s _1 cm -2 
within 6 years. 

Subject headings: gamma-rays: bursts -gravitational waves - cosmology: miscel- 
laneous 



1. INTRODUCTION 



The brightness distribution of cosmological sources is conventi onally used to constrain 



the luminosity function of the sourc es, their evolut i on in density (Peebles 



transient sources, their rate density (ISchmidt 



2001 



Sethi fc Bhargavi 



2001 



9931) and, fo r 



Totani 



1993). 



This method is applicable both to long-lived sources such as galaxies and to transient events 
such as supernovae and gamma-ray bursts (GRBs). Estimates are obtained by fitting the 
number - brightness distribution to models that include luminosity, source density and evo- 
lution effects. In the case of transient events an additional parameter is available - the event 
arrival times. 



The temporal distribution of transient astrophysical pop ulations of events h a s bee n 



described by the 'probability event horizon' (PEH) concept of 



Coward fc Burmanl (120051 ). 



This method establishes a temporal dependence by noting the occurrences of successively 
brighter events in a time series. By utilizing the fact that the rarest events will preferentially 
occur after the longest observational periods, it produces a data set with a unique statistical 
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signature. 

Here we show that a well-defined observation-time dependence is an intrinsic feature of 
the source distribution of events. Using Swift GRB data we demonstrate how this property 
can be used to constrain source distributions. We start by presenting an analytical derivation 
of the peak flux-observation time relation, P(T), for sources which are uniformly distributed 
in Euclidean space and then describe its extension to cosmological models ($2D- We derive a 
simple power-law relation for P(T) that is invariant to the luminosity distribution of events. 

We then utilize the PEH technique to show how P(T) data can be extracted from a 
distribution of peak fluxes (§3]). We show that the PEH method can be greatly improved by 
invoking time reversal invariance and the temporal cosmological principle: for time scales 
that are short compared to the age of the Universe, a distribution of independent events 
is invariant with respect to temporal direction and there is nothing special about the time 
when we switch on our detector. 

As a test, we apply the P(T) relation to Swift long-GRB data. We demonstrate that 
the technique can be used as a probe of the event rate density and luminosity distribution of 
the sources (Q, and as a tool to predict the likelihood of future high peak flux events (§H]). 

2. THE PEAK FLUX - OBSERVATION TIME RELATION 

In this paper we will define an event to be an astrophysical transient occurrence with 
a duration much less than the period of observation. Examples are GRBs and gravitational 
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wave burst sources such as coalescing compact binaries or core-collapse supernovae. 

Consider a distribution of events defined within a Euclidean space by an event rate 
density r and a luminosity functional (f>(L) (L min < L < L max ). The observed peak flux, or 
'brightness', distribution of events over an observation time T is a convolution of the radial 
distribution of the sources and their luminosity function. For peak fluxes (photons cm _2 s _1 ) 
between P and P + dP: 

&N(P)=AixT I (j)(L)dLr r 2 dr, (1) 

with r = (L/47TP) 1 / 2 . The total number of events observed in time T with a peak flux 
greater than P is given by : 



•imai fy/ L/A-kP 

N(>P) = TAnl <f>(L)dL r r 2 dr, (2) 



L m i n J 

where the average solid angle covered on the sky has been accounted for by Af2/47r. The 
upper limit in the integration over r is the maximum distance for which an event with 
luminosity L produces a peak flux P. 

For r and <j){L) independent of position, integrating over the radial distance yields: 



: We use here the luminosity function for GRB sources, 4>{L), which includes a normalization constant to 
ensure that it integrates to unity o ver the range of source lu minosities. This means that <fr(L) has units of 



inverse luminosity-see for example. 



Porciani fc Madaul POOII ). 
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N{>p) = Tr^W* p - 3/2 

3V47T J L 



<P(L)L 3/2 dL . 



(3) 



This is the familiar log iV-log P relation, N(> P) oc P 3 ^ 2 , a pow er law independent of the 



form of the luminosity function (IHorack. Emslie fc Meegan 



1994 ). 



To introduce the temporal distribution of events, we note that, as the events are in- 
dependent of each other, the individual events will follow a Poisson distribution in time. 
Therefore, the temporal separation between events will follow an exponential distribution, 
defined by a mean event rate R{r) = r (4/3)7rr 3 for events out to r. The probability for at 
least one event > P to occur in a volume bounded by r during an observation time T at 
constant probability e is given by: 



V(n > 1; R(r),T) = 1 - e R(r)T = e . 



(4) 



For this equation to remain satisfied with increasing observation time: 



N(> P) = R(r)T = |ln(l - e) 



(5) 



Equations (131) and (jSJ) for N(> P) combine to give the relation for the evolution of 
brightness function of observation time: 



P(T)- 



3V47r|ln(l-e) 



2/3 



(f(L)L 3/2 dL 



2/3 
j>2/3 



(6) 
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This relation shows that for a simple Euclidean geometry, a log P-log T distribution will 
have a slope of 2/3, independent of the form of the luminosity function. One can consider 
that changes in r$ create a horizontal offset in the log P-log T distribution, while changes 
in the integrated luminosity create a vertical offset. However, the slope is fixed by the 3-D 
Euclidean geometry. 

We can use the log P-log T relation to produce curves defining the probability, e, of 
obtaining some value of peak flux, P, within an observation-time, T, for a given r and <f>(L). 

For a cosmological distribution of sources, equation (jSJ) must be modified to allow for cos- 
mic evolution. A standard Friedman cosmology can be used to define a differential event rate, 



dR(z), in the redshift shell z to z+dz. The 
l umin osity distance dl(z) (see for example 



uminosity and flux will be rel ated through z by a 



Coward fc Burmanl (120051 ) or 



Porciani fc Madau 



(120011 )). In this case, solving equation (jSJ) numerically, with P = L/47rd^(z), will yield the 
cosmological log P-log T relation. 



3. AN ENHANCED PEH FILTER 



Coward fc Burman 



b utilize the time domain, we use the probability event horizon (PEH) filter of 
(120051 ) to produce P(T) time data. The PEH filter is a tool that exploits the temporal in- 
formation encoded in a time seri es of transien t even ts and works by recording successively 
brighter events in a time series. 



Howell at al. 



(120071 ) demonstrated that the unique statis- 



tical signature of events filtered in this way could be exploited to obtain rate estimates of 
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transient events. However, the significant probability of a bright event occurring early in an 
observational period meant that only a small fraction of data was used by the method. As 
a result, large uncertainties were obtained in the estimates. There are however, two ways in 
which the amount of usable data can be increased. 

Firstly, the temporal cosmological principle implies that the PEH signature of a transient 
population of events is independent of when a detector is switched on. Secondly, time reversal 
invariance allows the PEH filter to be applied to a data set in both temporal directions. Thus, 
a time series of events can be treated as a closed loop which can be interrogated in both 
directions. The observational period is now defined as the total length of the loop. The start 
time for the PEH analysis is now arbitrary so, without loss of generality, we can choose any 
start time. This allows the PEH filter to be applied in such a way that the brightest event 
can be set as the final event in a series. This ensures that the PEH filter is applied to the 
full data stream and the process can be repeated in each direction, increasing the quantity 
of PEH data. We refer to these techniques as 'from max' plus 'time reversal' (FMTR). We 
show below how FMTR increases the PEH sample and significantly improves the statistical 
resolution when applied to the Swift data. 



4. APPLICATION TO SWIFT DATA 

In this section, we will apply the log P-log T relation to a cosmological population of 
long GRBs. To account for the event rate and luminosity function, we will use estimates 



based on recent studies. Using the FMTR method, we will extract a time-dependent sample 
of GRB peak fluxes from the Swift data, and demonstrate how it can be constrained by a 
log P-log T fit. 

For our long-GRB peak flux sample we use data recorded by the Swifio satellite between 
2004 December and 2007 April. We consider only bursts with confirmed peak fluxes detected 
within the 15-150 keV band of the Burst Alert Telescope (BAT) and with T90 > 2s (a T90 
duration is the interval in which a signal contains 90% of its total observed counts). The 
total sample consists of 190 peak fluxes. 

Figure [TA displays the Swift peak flux distribution of long GRBs as a time series. It is 
apparent that as observation time increases, there is an greater probability of a bright event. 
By extracting successively brighter events as a function of observation time, the PEH filter 
samples events from the low probability tail of the distribution. The strong brightness-time 
dependence of these events creates a unique statistical signature which can be modeled by 
the log P - log T relation. 

In Table 1 we show the PEH filtered data. It is apparent how the time intervals between 
successive events increase with observation time. This is a result of a progressive sampling 
of the rarer events of the distribution. 

To apply the log P-log T relation to the filtered data we must first set up a model 
to account for the source rate evolu t ion a nd luminosity distribution. We use a model from 



the recent study of iGuetta fc Piran 



(120071 ). They employ a 'flat-A' cosmological model and 



2 This data can be obtained from the Swift website 



http: / /swift. gsfc.nasa.gov/ docs/swift/archive/grb_table.html 
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Hq = 65 km s 1 Mpc 1 for the Hubble parameter at the present epoch. For the isotropic 



luminosity function o 
based on the work of 



GRB p e ak lu minosities, 4>(L), they use a broken power law form 



Schmidt 



(120011 ). Assuming that the rate of GRBs traces the global 
star formation history of the Universe, they employ a number of different star formation 
rate models. For each, they determine best fitted values for the luminosity function and 
event rate density. For this stud y, we use their mode l (i) pa rameters, which are based on the 



SF2 star formation rate model of iPorciani & Madaul (120011 ) . These parameters include fitted 



values for the luminosity function and a local event rate density r = 0.1 Gpc 3 yr 1 . To 
account for the average solid angle covered on the sky by Swift, we use a value of AQ = 1.33 



(IBand 



20051 ). 



Figure [T]B shows a lo g P log T fit to the P EH filtered data (shown by squares), using 



the fitting parameters of 



Guetta Piranl (120071 ): i.e. there are no free parameters in the 



comparison between theory and observation. We define a 90% confidence band - shown by 
the shaded area - corresponding to the e = 95% (top) and e = 5% (bottom) probabilities 
of detecting at least one event within an observation time T. We see that the data is well 
constrained. The fit shows that by using only a small sample of the brightest events, it is 
possible to extract the geometrical signature of the source population and to test estimates 
of the luminosity function and rate density of events. 

The dashed lines of Figure [TJ3 show the 90% confidence band corresponding to the 
Euclidean model using equation (jSJ). We see that two of the first few events lie outside the 
Euclidean curves but are constrained by the cosmological model. These events, occurring 
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at early observation times, most likely result from sources at large cosmological distances. 
The very bright events at late observation times are more probable - it is apparent that the 
Euclidean and cosmological curves begin to converge in this regime. To account for non- 
uniformly distributed sources, the method could be refined to take into account the spatial 
distribution of potential host galaxies. 

To test the power-law dependence in equation (jSJ), we have performed least-squares 
fitting to the PEH filtered data using the Euclidean log P log T curve as a linear regression 
model. By setting the power as a single free parameter, we obtained a value of 0.67 ± 0.02, 
confirming the 2/3 slope derived in section [2j 

5. A PREDICTIVE APPLICATION OF THE LOG P - LOG T RELATION 

Figure [2] illustrates how the log P-log T relation can be used as a predictive tool. By 
mapping the temporal evo. ution of datect ionptobabiHty, the maximum brightness of fntnre 
events can be constrained (jCoward at al.ll2005l ). As in Fig. [TB, the shaded areas show the 
log P-log T detection confidence bands corresponding to different probability values (shown 
in the legend). The current Swift observation time (398 days) is shown by the vertical solid 
line. For predictive applications it is essential that the true temporal sequence of events is 
retained. Therefore, rather than using the FMTR technique, we apply the unmodified PEH 
filter from the time of the first event. Comparing with Fig. [TB, we see that the FMTR 
method has increased the sample by 220%, of which 80% is gained by incorporating time 
reversal. Using 100 Monte-Carlo simulations, we find a mean fractional increase in data of 



~ 35%, compared to the unmodified PEH filter. The larger than expected filter output using 
Swift data, implies that the FMTR method can be further optimally tuned. 

The plot shows that after 3 days of operation, there was a < 5% probability of detecting 
an event with peak flux equal to the first in the PEH sample, GRB 041123. The Poisson 
probability of detecting at least one event within a Gpc at this time is ~ 2 x 10 -6 . This 
implies that this event occurred at a considerable cosmological distance. The event GRB 
050525A is the brightest long GRB with a secure redshift, z = 0.606 - the Poisson probability 
of at least one event within this volume is 35%. The next brightest event, GRB 060017, is 
the most intense burst, in terms of peak flux, detected by Swift. 

As a demonstration of the predictive nature of the the log P - log T relation, Figure [2] 
shows that there is a 50% probability of obtaining an event with a peak flux greater than 
that of GRB 060017 within the next year and an 80% probability within 5 years. 

The curve predicts that there is a 50% (80%) probability of obtaining a burst with peak 
flux > 60 photons s -1 cm~ 2 within 2 (8) years. To determine the feasibility of this prediction, 
we consider again GRB 050525A, which had a peak flux of 42.3 photons s~ 1 cm~ 2 . Using this 
burst's redshift and converting to a peak luminosity, we find that an equivalent burst would 
have to occur within z ~ 0.5 to produce a peak flux of this level. The Poisson probability of 
at least one event at this distance within the next two years is ~ 90% (99%). If we consider 
a burst with a peak flux of 100 photons s _1 cm -2 , we find that an event of this peak flux is 
50% probable within 6 years. Such an event would correspond to a GRB 050525A-equivalent 



burst occurring within z ~ 0.4, for which the probability is 99%. 

The log P - log T technique naturally uses the brightest events of a data set. As the 
probability of obtaining a GRB afterglow increases with peak flux, the method can be used 
to predict the expected occurrence of events at low z. 

6. SUMMARY 

We have provided a clear demonstration of the log P-log T relation by applying the 
PEH filter to the Swift GRB peak flux distribution. A log P-log T model with no free 
parameters was fitted to filtered data confirming the power law P oc T 2 / 3 in the Euclidean 
limit; the power law is independent of the form of the luminosity distribution. 

The FMTR method significantly improves the PEH method, which was previously dis- 
advantaged by using only a small fraction of a data stream. We have shown that FMTR 
enables the PEH filter to use over 8% of the available data, making it a practical tool for 
cosmology. 

We have shown that the PEH technique can be used as a predictive tool. Comparing 
observation with prediction provides an additional means to test rate estimates and evaluate 
source parameters such as the limits of the luminosity distribution. 

In a future study, we intend to apply the FMTR method to both the Swift and BATSE 
GRB data. We will investigate the efficiency of the method in determining constraints on 



the rate density and limits of the luminosity function. 
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GRB 


Peak Flux 


Redshift 


Observation Time 




(photons s~ 1 cm -2 ) 




(days) 


060202 


0.5 




8 


060203 


0.6 




9 


060204B 


1.3 




10 


060206 


2.8 


4.045 


12 


060223B 


2.9 




29 


060306 


6.1 




42 


060418 


6.7 


1.49 


84 


060510A 


17.0 




106 


061121 


21.1 


1.314 


297 


050219B 


25.4 




506 


050525A 


42.3 


0.606 


602 


060117 


48.9 




839 


060111A 


1.72 




4 


060110 


1.9 




5 


060105 


7.5 




10 


060603 


27.6 


2.821 


227 



Table 1: Data extracted from the Swift peak flux distribution of long GRBs using the PEH 
filter. The observation time is determined by treating the time series as a closed loop and 
setting the last event in the series to be the brightest burst. The lower set of data was 
obtained by invoking time reversal. A log P - log T fit to this data is shown in Fig. [T] 
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Fig. 1. — Panel (A) shows the Swift peak flux distribution as a time series. It is evident that 
as observation time increases the probability of a bright event increases. Panel (B) uses a 
PEH filter to extract the geometrical sign ature of the GRB dist ribution (shown as squares). 



Assuming an event rate of 0.1 Gpc 3 yr 1 (IGuetta fc Piran 



20071 ) , the shaded area shows the 



cosmology dependent log P-log T model corresponding to a (5 - 95)% confidence band. The 
equivalent model for a Euclidean geometry is shown by the dashed curves - the two outliers, 
which have no associated redshifts, probably result from distant cosmological events. 
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Fig. 2. — The log P-log T relation used as a predictive tool. The successive maximum peak 
fluxes detected by Swift since the start of operation are shown by squares. The shaded areas 
show detection confidence bands corresponding to different probability values (shown in the 
legend). The current Swift observation time (398 days) is shown by the solid line. The plot 
shows that after 3 days of operation, there was a < 5% probability of detecting an event 
with peak flux equal to the first in the PEH sample, GRB 041123. We see that within the 
next year, there is a 50% chance of obtaining a peak flux greater than for GRB 060017 - the 
most intense burst (in peak flux) detected by Swift - and the same probability of obtaining 
a burst with peak flux > 100 photons s~ 1 cm~ 2 within 6 years. 



